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Abstract 

We investigate the dynamics of two interacting electrons in coupled quantum 
dots driven by an AC field. We find that the two electrons can be trapped in 
one of the dots by the AC field, in spite of the strong Coulomb repulsion. In 
particular, we find that the interaction may enhance the localization effect. 
We also demonstrate the field excitation procedure to generate the maximally 
entangled Bell states. The generation time is determined by both analytic and 
numerical solutions of the time dependent Schrodinger equation. 
PACS numbers: 03.65.Ud, 78.67.Hc, 73.23.-b 
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I. INTRODUCTION 

Quantum-state engineering via optical or electrical manipulation over the coherent dy- 
namics of suitable quantum-mechanical systems has become a fascinating prospect of modern 



physics. A very intriguing result in the study of a single particle in a double-trap system 
exposed to a time-dependent external field is the dynamic localization phenomenon 
i.e., for appropriate field parameters, a localized wave packet remains dynamically localized 
during the subsequent time evolution. This driving induced trapping has been theoretically 
and experimentally studied in many physical and chemical systems j|,|5||J. 

When two or more interacting particles are present, apart from the highly nontrivial 
question of whether the strong many-body interaction can be overcome for the particles 
to create and preserve localization, the possibility of entanglement of the many-body wave 
functions arises. Entanglement is an essential ingredient in any scheme of quantum infor- 
mation processing. Recently, quantum dot realizations of the entanglement have received 
increasing attention. Various schemes based on electron spins and electron-hole pairs have 
been proposed to implement quantum computer hardware architectures [ff||8|,P, [10| , |TT| . Al- 



though there have been some numerical studies on interacting electron systems driven by 
an AC field [|l^,|T^,[Ll]] , there is still little theoretical understanding of the observed effects 
beyond the phenomenology level at the present time. 

In this paper, we address the dynamic localization and quantum entanglement of two 
interacting electrons in a double quantum dot system driven by an AC field (see Fig. 1). 
We show that with certain choice of parameters, in contrast to direct intuition, the Coulomb 
repulsion may enhance the localization due to the level crossing associated with different 
symmetries of a dynamic parity operation (discussed in more detail later). We also show 
that the maximally entangled Bell state can be prepared and maintained with a pulse of 
an AC field. Our study provides useful information for the future exploitation of coherent 
control of two-electron states in quantum dots. 

In Sec. II we present the Hamiltonian for two interacting electrons in a double-dot system. 
A simplified spin-1 representation of the Hamiltonian is given in Sec. III. In Sec. IV we 
discuss the phenomenon of Coulomb-enhanced localization. The field excitation procedure 
to generate the maximally entangled Bell state is showed in Sec. V. A summary is given in 
Sec. VI. 



II. THEORETICAL MODEL 



The Hamiltonian which we use to describe the dynamics of two interacting electrons in 
a coupled quantum dot driven by electric field is 

H(t)= £ ^p^+VbGn-i^l), (la) 

i=l,2 

h(r, p, t) = ^ - ezF(t) + V t (r) +V l (r), (lb) 

Vc(\t 1 -t 2 \)= e% (lc) 
k\x\ - r 2 | 

where VcQri — r 2 |) is the Coulomb interaction and h the single-particle Hamiltonian. The 
time-dependent electric field F(t) is applied along the inter-dot axis. The dielectric constant 
K and the effective mass m* are material parameters. The potential Vt in h describes the in- 
plane confinement, whereas Vi models the longitudinal double-well structure. The transverse 
coupling of the dots is modeled by a harmonic potential 

V t (x,y) = ^(x> + y*). (2) 

It has been shown in experiments with electrically gated quantum dots in a two-dimensional 
electron system that the electronic spectrum is well described by a simple harmonic oscillator 
1 375 1 - In describing the confinement V\ along the inter-dot axis, we use a (locally harmonic) 



double well potential of the form 

V t = ^4^-ar, (3) 

which, in the limit of large inter-dot distance, separates into two harmonic wells (one for each 
dot) of frequency Although in principle a square- well potential would be a more accurate 
description of the real potential than the harmonic double well, there is no qualitative 
difference between the results presented below obtained with harmonic potential and the 
corresponding results using square-well potential. 



We assume the transverse confinement is strong enough. Thus, the (x, y) degrees of 
freedom are frozen in the dynamics and the two-electron wave function can be written as 
|^(ri, r 2 , £)) = \4>(xi,yi))\(j)(x 2 ,y 2 ))\§( Z i, z 2 ,t)), where \<f>(x,y)) is the transverse ground 
state. Considering the fact that the external electric field is applied along the longitudinal 
direction, the approximation of frozen in-plane motion is reasonable and has been employed 
in previous work [13|]. After integrating over the (x,y) degrees of freedom, the time depen- 
dent Schrodinger equation becomes 

.d\<f>(z 1 ,z 2 ,t)) . H 2 ( d 2 d 2 \ rr , , rr , . 

'* m = [ ~2^~* \d? 1 + dz 2 ) + V{zi) + V{Z2) (4) 

+V c (\z t - z 2 \) - eF(t)( Zl + z 2 )]|$0i, Z2,t)), 

where Vc is the effective one-dimensional Coulomb interaction 

V c {\ Zl - z 2 \) = / dx 1 dx 2 dy 1 dy 2 , , . (5) 

J k\ti — r 2 \ 

We use the effective mass m* and dielectric constant k of GaAs. The other parameters are 
chosen as huii = 20meV and a = 20nm. 

The time evolution of the two electrons is obtained by numerically solving the Schrodinger 
equation (4) using an extension of the Crank-Nicholson method [|T(J to two spatial dimen- 
sions. The initial state used in this investigation is the field-free ground state, which is 
obtained by propagating a trial two-particle wave packet in the imaginary time domain and 
shown in Fig. 2(a). As illustrated by the spatial symmetry under the particle exchange, the 
ground state is a singlet, as expected. Since Eq. (4) contains no mixing between the singlet 
and triplet sub-spaces, the dynamics of the system will always be confined to the singlet 
sub-space. 

III. SPIN-1 REPRESENTATION OF THE HAMILTON! AN 

To understand the underlying physics behind the numerical results we employ the 
Hund-Mulliken (HM) approximation by introducing the orthonormalized one-particle wave 
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functions |<&±o) = (\f±a) — g\ip Ta ))/ ~ 25g + g 2 , where \<p± a ) are the single-particle 
ground states for the right and left dots, S = (</? +a | </?_„) denotes the overlap integral, 
and g — (1 — y/l — S 2 ) / S. Using |$± a ), we construct three singlet basis functions with 
respect to which we diagonalize the two-electron Hamiltonian: Two states with double oc- 
cupation in each dot, \RR) = \<S> +a (zi)) \$ +a (z 2 )), \LL) = \^^ a (z 1 ))\^^ a (z 2 )), and one state 
with single occupation in each dot, \LR) = (l/ v / 2)[|$+a(^i)) \®-a(z 2 )} + \$ +a (z 2 )) |$_ a (2i))]. 
Calculating the matrix elements of the Hamiltonian in this orthonormal basis, we obtain 

H HM {t) = uJ 2 z + 2wJ x - n(t)J t , (6) 

where we have dropped a constant energy term that makes no contribution to the dynamics. 
In Eq. (6) u = {LL\Vc\LL} — (LR\Vc\LR) is the difference between the intradot and interdot 
Coulomb interaction, w = (&± a \h(z =F a)|$ Ta ) + (LR\V C \LL) denotes the single-particle 
tunneling induced by dot-dot coupling and the Coulomb interaction, //(£) = 2eF(t)a(l — 
g 2 )/(l — 2Sg + g 2 ) describes the electron-field coupling, and J x and J z are the x- and z- 
components of the spin-1 operator. The localized two-particle state \LL) is equivalent to 
the eigenstate | j z — 1) of J z and \RR) to the state | j z = —1), while the delocalized state 
\LR) is identical to the state | j z = 0). According to the expression for u and w, using the 
material parameters given below Eq. (5) we obtain u = 12meV and w = — 0.4meV. 

In the absence of an external driving, the eigenenergies of Hum are E\ = (u — 
\/u 2 + 16u> 2 )/2, E 2 = u, and E 3 = (u + \/u 2 + 16w 2 )/2, and the corresponding eigenstates 
are 

|v?i> = \LL) - E 3 /(V2w)\LR) + \RR), (7a) 
M = (\RR) - \LL))/V2, (7b) 

|^ 3 > = \LL) - E 1 /(V2w)\LR) + \RR). (7c) 

The symmetric ground state |</?i) is dominated by the delocalized state \LR) due to the strong 
Coulomb repulsion, whereas the other two eigenstates are nearly degenerate and dominated 
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by the two localized states \LL) and \RR). The superposition of the two localized states \LL) 
and \RR) implies that the spatial wave functions of the two electrons have been entangled 
and correlated, in the usual sense that they cannot be factorized into single-particle states. 
The nonlinear term in H HM can be exploited to generate entangled states. The ground state 
of Eq. (7) is plotted in Fig. 2(b). Clearly the HM approximation describes the ground state 
very well when compared with the exact numerical solution [Fig. 2(a)]. 

IV. COULOMB-ENHANCED DYNAMIC LOCALIZATION 

To investigate dynamic localization, one must start with a localized wave packet. This 
can be realized from the unperturbed ground state by two separate methods: One is to 
suddenly switch on a DC field with the strength satisfying the resonance condition /i = u. 
At time t ps n/(y/2w), the delocalized ground state l^) evolves into a localized state \RR) 
with two electrons occupying the same right dot. The alternative method is to adiabatically 
switch on a constant electric field. Then as shown in Fig. 2(c), the ground state configuration 
develops a series of Coulomb stairs as a function of the field amplitude, corresponding to a 
series of avoided crossing in the energy spectrum [see Fig. 2(d)]. 

After preparing a localized state, say \RR), by adiabatic evolution in the presence of a DC 
field, the DC field is suddenly swiched off and an AC field of the form F 1 sm(ut) is swiched 
on. We show now that the localization can be dynamically maintained by the AC field 
even when the field strength is small compared to the Coulomb repulsion between the two 
electrons. Time periodicity of the Hamiltonian enables us to describe the dynamics within 
the Floquet formalism. In addition, since the Hamiltonian is invariant under the combined 
dynamic parity operation z — > —z; t — > t + tt/cu, each Floquet state is either odd or even. 
Quasienergies of different parity may cross, otherwise an avoided crossing may occur. Figures 
3(a)-(b) show the quasienergy spectra of HHuif) as a function of iq with the presence and 
absence of the Coulomb interaction, respectively. Two prominent features can be identified 
from the comparison of these two cases: (i) the strong Coulomb interaction removes the 



6 



level crossings among three two-particle states and thus induces avoided crossings; (ii) in 
the weak field regime, there occurs a crossing between the quasienergies e 2 and e 3 , which 
develop from the unperturbed eigenenergies E 2 and E 3 . To illustrate the effect of this 
crossing on the system's dynamics, we begin with the initial state \RR) and follow the time 
evolution of the probability PRu(t) for finding the two electrons in the right dot. The result 
is shown in Fig. 3(c). It is clear that Prr always remains near 1 as if the two electrons were 
frozen in the same right dot. This dynamic localization seems counterintuitive at first as the 
Coulomb repulsion is very strong compared to the AC field, preventing the two electrons 
occupying the same dot. The localization shown in Fig. 3(c) has no correspondence in a 
non-interacting case, since no level crossings occur in the weak field regime if the Coulomb 
interaction is absent. Therefore, the present localization effect results from an interplay of 
the Coulomb interaction and an AC field. Furthermore, we find that even if the Coulomb 
interaction is strong enough, the dynamic localization can still occur at the crossing of the 
quasienergies e 2 and e 3 . Slightly away from the crossing point, the life time for finding 
two electrons in the right dot is very long. This makes the experimental realization of the 
dynamic localization more realistic. The ratio 2/ii/uj at the first crossing is about 2.4, a 
root of the zero-order Bessel function, suggesting that the two-electron localization can be 
approximated by a two-state model composed of \ip 2 ) and | v?3 ) . 

Surprisingly, we find that in the strong field regime, the presence of Coulomb repulsion 
may help to enhance the dynamic localization when compared with the non-interacting case. 
From Eq. (7) we see that due to the strong Coulomb interaction, \RR) ~ (l/^/2)(\tp 2 ) + \(p 3 )). 
Thus the initial localized state can be approximated by a superposition of the degenerate 
Floquet states developed from the nearly-degenerate states \(p 2 ) and |</? 3 ), which leads to 
a complete suppression of tunneling at the crossing of the quasienergies e 2 and £3 even 
in the weak field limit. In the absence of the Coulomb interaction, the initial localized 
state \RR) is a superposition of all three eigenstates, suggesting dynamic localization at the 
crossing among three quasienergies. However, the fundamental difference lies in the fact that 
when the strong interaction is present, the tunneling coupling is (LL\Hhm\RR) — 4:W 2 /u, 
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whereas in the absence of the Coulomb interaction the tunneling coupling is 2w. Thus 
the coupling between the two localized states greatly decreases in the presence of strong 
Coulomb interaction, which leads to the enhancement of localization effects. 

We emphasize that the localization only occurs at the crossing of the quasienergies de- 
veloped from the high-energy eigenstates |</? 2 ) and |</? 3 ), in which case the delocalized ground 
state \ipi) is a dark state during time evolution. For comparison, Figure 3(d) shows Pim(t) 
with the field parameters corresponding to the crossing of the quasienergies E\ and £3. It 
reveals that dynamic localization disppears and the two electrons start to oscillate between 
two dots. This fact highlights the essential difference between a single-particle system and an 
interacting two-particle system. The level crossing of S\ and e 3 induces strong participation 
of the Floquet state evolved from the unperturbed ground state. As its largest component is 
the delocalized two-particle state \LR), dynamic localization will be completely destructed, 
as shown in Fig. 3(d). 

V. BELL STATE GENERATION 

In this section we investigate the field excitation procedure to obtain a maximally en- 
tangled Bell state of the form \^Beii) — (\RR) + e^|LL)) with an arbitrary phase angle (p. 
Figure 4(a) plots the Floquet quasienergies of H HM (t) as a function of the driving frequency 
uj. It shows that when uj = u, an avoided crossing is formed between the quasienergies e± 
and £2- Choosing the field parameters at this avoided crossing and starting with the un- 
perturbed ground state, we present in Fig. 4(b) the probabilities PLFi(t) for finding the two 
electrons in different dots, Pll in the left dot, and Prr in the right dot. It shows that the 
two electrons oscillate between the delocalized state \LR) and two localized states \LL) and 
\RR). The occupations of two localized states are always the same and the oscillations are 
always in-phase. Figure 4(c) shows the probability pseii for finding the maximally entangled 
Bell state with = it. One can see that pseii varies with time, reaching a maximum value 
of 1 when P LR = (the two electrons are maximally entangled with p BeU = 1). A direct 
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numerical solution of Eq. (4) gives the same prediction. 

Once the two electrons are in the maximally entangled Bell state, they can remain 
maximally entangled by suddenly turning off the AC field. We show this effect in Figure 
5. Figure 5(a) plots the time evolution of the occupations of three two-particle states, and 
Figure 5(b) the probability p BeU with = n (solid line). It is clear from Fig. 5(a) that a 
pulse of an AC field induces the two electrons to stay in the same dot, while each of them 
occupies either of the dots with the same probability. As shown in Fig. 5(b) (solid line), the 
two electrons remain maximally entangled after the electric field is turned off. Therefore, 
the maximally entangled Bell state can be created and maintained by applying a pulse of a 
resonant AC field. 

We derive an approximate analytical solution to highlight the physical aspect of the 
Bell-state generation procedure. Taking into account symmetric properties of the three 
unperturbed eigenstates of the system, the dynamics shown in Fig. 5 is determined by the 
one-photon transition between the ground state \<pi) and the first excited state \<p 2 ), whereas 
the transition from \ifi) to \ip 3 ) is prohibited due to their identical symmetry. In this case, 
we can approximate the Hamiltonian H HM (t) in a Hilbert space spanned by the states |</?i) 
and \ip 2 ) 



H HM(t) 



1 E 1 -V2fi(t)/X ^ 
y -V2fi(t)/X E 2 



(8) 



where X = ^Aw 2 + E~$/y/2w. In the case of one-photon resonance E 2 — E 1 = to and after 
applying the rotating-wave approximation we obtain the expression for the probability of 
finding the maximally entangled Bell state 

PBeii (t) = ^ (1 - cos 4>) sin 2 (wfj.it/ u), (9) 

in the interaction representation and weak coupling limit u ^> w. We can see from Eq. 
(9) that the quantum state of the system at time r = Tru/(2w/ii) corresponds to a = n 
maximally entangled Bell state (\RR) — \LL))/\^2. 



The result of Eq. (9) is shown in Fig. 5(c) (dotted line). Clearly, in comparison with 
the exact numerical solution (solid line), our two-state approximation describes the system 
evolution very well. Therefore we arrive at the conclusion that a selective pulse of an AC field 
with duration r = nu/(2wfii) can be used to create and maintain a maximally entangled 
Bell state in the system of two electrons in two coupled quantum dots. 

We notice from Eq. (9) that Bell-state generation time is significantly shortened by 
increasing the amplitude of the AC field. This is important because shorter Bell-state 
generation time is fundamental to the experimental observation of such maximally entangled 
state, which is impeded by inevitable decoherence occurred in the realistic double quantum 
dot system. The decoherence is the most problematic issue pertaining to most quantum 
computing processing. In the present entangled state proposal, the decohering time depends 
partly on the fluctuation of the single particle energy caused by the modification of the 
confining potential due to phononic excitations. There is also a quantum electrodynamic 
contribution because of coupling to the vacuum modes. In addition, impurity scattering and 
phonon emission also have contributions to the decoherence. However, in principle, their 
effects can be minimized by more precise fabrication technology and by cooling the system. 

VI. SUMMARY 

In summary, we have shown how the dynamic localization and entanglement of two in- 
teracting electrons in a coupled quantum dot system can be generated and maintained using 
an AC field. The Bell-state generation time has been calculated by an analytical approach. 
Exact numerical calculation confirms these two kinds of Coulomb-involved phenomena. We 
hope the present study will shed some light on the future development of the coherent control 
of electrons in quantum dot systems. 
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Figure Captions 
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FIG. 1. Sketch of two interacting electrons in a coupled quantum dot driven by electric 
fields. 

FIG. 2. (a) Two-electron ground state versus the longitudinal coordinates of the two 
electrons along diagonal direction z\ = —z 2 = z, obtained via numerically integrating Eq. 
(4) with no AC field; (b) Two-electron ground state obtained analytically via the Hund- 
Mulliken approximation; (c) Electron-number distribution of the right dot in the ground 
state as a function of DC electric field; (d) Low-energy spectrum versus DC field. 

FIG. 3. (a) Quasienergy spectrum for two interacting electrons for uj = 4.01meV; (b) 
Quasienergy spectrum with no Coulomb interaction; (c) Probability Prr that two electrons 
occupy the same right dot for an AC field of F 1 = 1.078kV/cm, which produces an ex- 
act crossing between two quasienergies that developed from nearly-degenerate two-particle 
eigenenergies; (d) Prr for F 1 = 0.76kV/cm and uj = 6.4meV, corresponding to a crossing 
between two quasienergies, one of which originates from unperturbed ground-state. 

FIG. 4. (a) Floquet spectrum versus the driving frequency of an AC field with a strength 
Fi = lkV/cm; (b) Probabilities that the two electrons are in different combinations of two 
quantum dots, for a driving frequency uj = 12meV, corresponding to the avoided crossing 
shown in (a); (c) Probability pseii as a function of time. 

FIG. 5. (a) Probabilities that the two electrons are in different combinations of two 
quantum dots, with the AC field switched off after t = 8.1ps; (b) Probability p Be ii as a 
function of time. The numerical and analytic results are shown with the solid and dotted 
lines, respectively. 
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